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We consider continuum dielectric models as minimal models to understand the effect of a surround- 
ing globular protein and solvent on the quantum dynamics of electronic excitations in a biological 
chromophore. For these models we derive expressions for the frequency dependent spectral den- 
sity which describes the coupling of the electronic levels in the chromophore to its environment. 
The magnitude and frequency dependence of the spectral density determines whether or not the 
quantum dynamics is coherent or incoherent, and thus whether on not one can observe quantum 
interference effects such as Rabi oscillations. We find the contributions to the spectral density from 
each component of the chromophore environment; the bulk solvent, protein, and water bound to 
the protein. The relative importance of each component to the quantum dynamics of the chro- 
mophore is determined by the time scale on which one is considering the dynamics. Our results 
provide a natural explanation and model for the different time scales observed in the spectral density 
^ I extracted from the solvation dynamics probed by ultra-fast laser spectroscopy techniques such as 

^ ■ the dynamic Stokes shift and three pulse photon echo spectroscopy. Our results are used to define 

under what conditions the dynamics of the excited chromophore is dominated by the surrounding 
1/^ ' protein and when it is dominated by dielectric fiuctuations in the solvent. We show that even when 

^SJ , the chromophore is shielded from the solvent by the protein ultra-fast solvation can be dominated 

by the solvent. Hence, we suggest that the ultra- fast solvation recently seen in some biological 
' chromophores should not necessarily be assigned to ultra-fast protein dynamics. The magnitudes of 

^ , the spectral density that we estimate from our continuum models and extracted from experiment 

'/^ ■ suggest that most quantum dynamics of electronic excitations is incoherent. A possible exception 
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I. INTRODUCTION 



The functionality of many proteins is associated with a small subsystem or active site such as a heme group, a couple 
of amino acids involved in proton transfer, or a co-factor such as an optically active molecule (chromophore). There 
are a diverse range of optically active molecules that have an important biological function Jj. Examples include 
^ ■ retinal (involved in vision), green fluorescent protein and porphyrins (photosynthesis). For these chromophores, the 
' protein acts as a transducer which converts optical excitation of the chromophore into a change such as an electrical 
^ signal or conformational change that in turn brings about the desired biological function. Many of these transducers 

operate with speeds, specificities, and efficiencies which nanotechnologists are striving to mimic 
^ , The dynamics of a protein involves thousands of degrees of freedom and at room temperature can be described by 
classical mechanics and modelled using molecular dynamics methods. In contrast, the functional subsystem involves 
only a few quantum states and their dynamics must be described quantum mechanically. This has led to considerable 
effort at developing hybrid QM/MM (quantum mechanical-molecular mechanical) methods [1, In niost cases 
the change in quantum state associated with the functional event (e.g., transition to an excited electronic state) is 
associated with a change in the electric dipole moment of the subsystem. Since the protein contains polar residues 
and is surrounded by a highly polar solvent (water) J,, ^ there is a strong interaction between the functional 
subsystem and its environment. Consequently, the environment can have a significant effect on the quantum dynamics 
of the subsystem. Indeed chromophores such as retinal, Photoactive Yellow ProteinQ, Green Flourescent Protein 
exhibit distinctly different dynamics in solution, in the gas phase, and in the protein environment. [^,[l0| For example, 
the speed, efficiency, and selectivity with which excited retinal undergoes a conformational change are all significantly 
less in water than in the protein environment. 0, [TTI . [T^ This interplay between quantum and classical dynamics 
raises a number of questions of fundamental interest. On what length and time scales does the crossover from quantum 
to classical behaviour occur? When are quantum mechanical effects such as coherence (i.e., superposition states), 
entanglement, tunneling, or interference necessary for biological function? [l^, [ij] What aspects and details of the 
structure and dynamic properties of the protein are crucial to biological function? 
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FIG. 1: The chromophore, protein and bound water in photoactive yellow protein (PYP). The isolated spheres represent the 
bound water, the chromophore is shown by its van de Waals surface, and the protein by a cartoon representation. Observe that 
the protein surrounding the chromophore reduces the contact of the chromophore with the surrounding bulk water. Generated 
from the Protein Database 3PYP.pdb[2||. 



A. Biomolecular chromophores 

Most chromophores are large conjugated organic molecules which are surrounded by a protein which in turn is 
surrounded by a solvent. Figure [T] shows the photoactive yellow protein (PYP), including the chromophore and 
the so-called "bound water" molecules which reside with comparatively long lifetimes on the surface of the protein 
p^ . Most chromophores have large dipole moments which change significantly upon optical excitation, leading to 
significant relaxation of the polarisable environment. Combined chromophore- protein-solvent systems exhibit a broad 
range of time, length and energy scales (see Figure [5]). Typical values of different time scales are shown in Table 1 in 
the Appendix. 

In this paper, we focus on minimal model Hamiltonians, since we are seeking to understand general qualitative 
features, identify crucial parameters for understanding qualitative changes in behaviour. We specifically consider 
chromophores which can be described as two level systems (TLS), e.g., we only need consider the ground and first 
excited state. The models proposed here can also be extended to include internal nuclear dynamics of the chromophore, 
such as conformational change 

Questions of quantum coherence and the role of the environment are particularly pertinent and controversial in 
photosynthetic systems [l^, [13, [3 [3| • It has sometimes been claimed that the excitons within the light harvesting 
rings are quantum mechanically coherent over some or all of the chromophores (sometimes as many as 32) within the 
ring. It has also been suggested that such coherence is important for optimum performance of the system [20l. 
On the other hand, inter-ring transfer of excitons is incoherent which ensures the desirable feature of irreversible 
transport of energy towards the reaction centre. Recently, the question of delocalisation of excitons over several base 
pairs in DNA has be studied, motivated by a desire to understand UV damage of DNA.[23| 

B. Quantum dynamics, decoherence, and the spin-boson model 

Understanding, the dynamics of a quantum system which is strongly coupled to its environment is a challenging 
theoretical problem that has attracted considerable attention over the past few decades Many consider that 

decoherence is the key to resolving the quantum measurement problem, [ssl [sgI . [stI These issues are receiving renewed 
interest because decoherence is detrimental to quantum information processing . l38l [39j Substantial progress has been 
made by considering the simplest possible models such as the spin boson model jSJ] which describes a two-level system 
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FIG. 2: Schematic representation of the time scales of various processes in proteins and solutions. ET stands for electron 
transfer, PS RC for photosynthetic reaction centre. Specific numbers and references can be found in the Appendix. 



TABLE I: Comparison of the matrix element A which couples two quantum states for various processes in proteins with the 
solvation rates due to the interaction of the quantum system with different parts of its environment. The quantum dynamics of 
the process will be determined largely by the part of the environment which undergoes solvation relaxation at a rate comparable 
to A. LHI and LHII refer to light harvesting complexes I and II in photosynthetic purple bacteria. 



Process 


A energy (meV) 


Ref. 


Forster coupfing between chromophores in FRET spectroscopy 


0.2-2 




Interring Forster coupling between chromophores in LHI and LHII 


0.3 


[20] 


Intraring Forster coupling between two chlorophyl molecules in LHI 


50-100 


[20] 


Forster coupling between infrared amide modes in proteins 


0.1-1 


[25] 


Electron transfer in photosynthetic reaction centre (PRC) 


1-10 


[26] 


Electron transfer in PRC radical cation 


10-100 


[22 


Electron transfer in proteins 


10"'' - 10-2 


[28] 


Proton transfer 


0.05 


[29] 


Level crossing for non-radiative decay 


40 


[30] 


Solvation rate due to bulk water 


10 


m 


Solvation rate due to bound water 


0.1-1 


[32] 


Solvation rate due to protein 


0.004-0.04 


[33] 



(the "spin")which is coupled linearly to an infinite "bath" of harmonic oscillators. 

For biomolecular systems, the spin-boson model has previously been applied to electron transfer. [J |43, El We have 
recently shown the relevance of the spin-boson model to understanding the effect of the environment on Forster reso- 
nant energy transfer between two chromophores Q. Of particular interest is the case where two molecules are coupled 
by Resonance Energy Transfer (RET), such as rings of chlorophyll molecules in photosynthesis and in Fluorescent 
Resonance Energy Transfer spectroscopy (FRET). Here, an excitation in one chromophore may be transferred to a 
nearby chromophore by the Coulomb interaction, typically dipole-dipole interactions. A coupled system of molecules 
such as this may be mapped to the spin-boson model [11, where the two quantum states refer to the location of the 
excitation, e is the difference in the two chromophore's excited energy levels, and J(w) describes the coupling of the 
excitation to the environments surrounding each molecule. We have previously shown [6] that the appropriate spectral 
density is simply the sum of the spectral density of each individual chromophore-protein complex. The magnitude of 
the spectral density then determines whether the transfer is coherent (oscillatory) or incoherent (one-way). There are 
several definite experimental signatures of the coherent interaction of a pair of chromophores. These include (Davydov) 
splitting of energy levels super- and sub-radiance (i.e., increase and reduction of the radiative lifetime ^22.. ,43j) 
and changes in fluorescence anisotropy 4J| . Both coherence (within a ring) and incoherence (between rings) may play 
potentially important functional roles in light harvesting complexes. 
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The Hamiltonian. For the oscillators this can be written as 



where (3 is an index denoting a particular oscillator with mass rrip, frequency LOp, momentum pfj, and position qp. 
We introduce second quantised operators and at where ap = \/— 11^(9/3 + "ipp/in^p^p)): a-nd satisfy the boson 



/3,7 



commutation relations [a^, a]y] = 

The Hamiltonian of the whole system is 

H = Ho + H2 (2) 

where 



In terms of the second quantised operators this can be written 

H = iecT^ + Aa^ + ^ fiuuiia^pap + (J^'^ C!,3{al + ap), (4) 

where and Uz are Pauli spin matrices, the C/3 describes the coupling of the system to each bath mode f3, e is the 
separation of system energy levels and A is the tunneling matrix element coupling the two states. Another model 
is the spin-bath model [i^, where the system of interest is coupled to specific localised states of the environment, 
themselves treated as two level systems. 

Spectral density. For the spin-boson model, the quantum dynamics of the two-level system (TLS) is completely 
determined by a single function, the spectral density, [sj] which is defined by 

47r 

It describes how strongly the oscillators with frequency near uj are coupled to the two-level system. An important 
quantity is the reorganisation energy defined by, 

Ej^^Tdu;^. (6) 



Many systems are described by ohmic dissipation, for which J{lu) ~ fiaoj, below some cutoff frequency, Wc, related to 
the relaxation rate of the environment, and if A ^ fiojc then at frequencies higher than this cutoff the coupling to 
the bath of oscillators can be neglected. The main purpose of this paper is to derive physically realistic expressions 
for this spectral density that are relevant to biological chromophores interacting with their environment. 

Known results. If e, A <C fiojc, for ohmic dissipation a is a critical parameter for determining the qualitative 
properties of the quantum dynamics [H, |4^. At zero temperature, for a < ^ the state of the TLS, exhibits damped 
Rabi oscillations, a signature of quantum coherence and interference. This can be described by considering the time 
dependence of the probability that the system is in one of the two levels, which can be related to the expectation value 
{az{t)). For i < a < 1, the system exhibits incoherent relaxation (exponential decay of {(Tz{t)), and for a > 1 the 
system is localised in its initial state - an example of the quantum Zeno effect (47j . A non-zero temperature reduces 
the range of a over which coherent oscillations can occur (see Fig. 21.2 in [33|). 

If A > hcjc then the results of Refs. [s^l and [i^ do not apply [4^. The Ijath responds slower than the relevant 
timescale for the dynamics of the TLS. Consequently, in order to destroy coherent oscillations, the bath must couple 
more strongly to the two- level system than for the case A <C hujc- System dynamics has been studied with quantum 
Monte Carlo simulations [49| when e = 0. Coherent oscillations in {az{t)) may be present for a > 1. Fig . 13 in Ref. 
(49[ shows that when A = Qfiujc then coherent oscillations can exist even for a = 30. Fig. 8 of Ref. [SO] shows that 
for A = TujJc, numerical renormalisation group calculations predict that coherent oscillations exist for a < 1.5. Fig. 3 
of Ref. shows that a renormalisation flow equation approach predicts that when A/hcUc increases from very small 
values to 0.3 that the allowed parameter range for coherent oscillations increases to a < 1.5. The coherent-incoherent 
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transition is associated with the delocalisation-localisation transition that has been studied widely in chemistry in 
systems such as the Creutz-Taube ion and the special pair in photosynthetic systems. j27i] 

The quantum dynamics of the spin boson model ([4]) for a general spectral density J(w) will be largely determined 
by the magnitude and frequency dependence of J{lu) for uj ^ A. For example, when the bath is weakly coupled (i.e., 
J(A) ^ A) to an unbiased (e = 0) two level system coherent oscillations exist, and the relevant decoherence rate 
given by Fermi's Golden Rule is[33| 

J_.j(A)„a(_^), (7) 



C. The chromophore environment: protein, bound water, and bulk water 

The structures and dynamics associated with the interaction of proteins with water is extremely rich and a challenge 
to model and to understand. 0, H, [13, [53| . One can classify the water molecules associated with proteins into several 
categories, (i) Water which is distant from the protein and has the same properties as bulk water, (ii) Water at the 
surface of the protein molecule. The first layer of molecules is referred to as the first hydration or solvation layer. 
These molecules are weakly bound to the charged residues found at the protein surface, (iii) Water buried inside the 
protein and which often binds to specific sites in the protein via multiple hydrogen bonds. The water inside and at 
the surface of the protein can exchange with the bulk water. 

Advances in experimental probes such as neutron scatteringfs^l, nuclear magnetic resonance (NMR)[8'], femtosecond 
laser spectroscopy, [5^] and dielectric dispersion [53| has allowed a quantitative description of the properties of the 
water molecules associated with specific parts of the solvated protein. Key quantities that can be determined include, 
(a) the occupancy (i.e., the probability that a water molecule will be found at the site), (b) the residence time (the 
timescale for exchange of the water molecule with the surrounding bulk water), and (c) the "order parameter" which is 
a measure of the rotational freedom of the water molecule at the site. NMR measurements suggest that the molecules 
at the surface exchange with the bulk water on timescales ranging from 10 psec to 1 nsec[8|. In contrast, buried 
molecules exchange with the solvent on timescales of the order of 1 ns to 1 /zs [sj . 

The term "biological water" has been used to describe water in proximity to a biological macromolecule [55| . 
Dielectric relaxation is significantly different in biological water [s^j . Whereas in bulk water the dominant dielectric 
relaxation time is 8.3 ps, for bound water this can be 2 to 4 orders of magnitude larger. Dielectric spectroscopy 
measurements of proteins in aqueous solutions found four dielectric relaxation times. [5J| For myoglobin these times 
were attributed to (1) reorientation of bulk water (8 ps), (2) relaxation of water associated with the protein (10 ps 
and 150 ps), and (3) reorientation of the whole protein molecule (15 nsec). 

Given the structural, chemical, and dynamical complexity of these environments we briefly discuss the limitations 
of some of the underlying simplications and approximations we assume in our models. Although, these simplifications 
may lead to quantitative differences between the predictions of our models and real systems we do not anticipate 
qualitative differences. 

Spherical symmetry. We assume that the chromophore is located at the centre of a spherical cavity inside a spherical 
protein. Similar assumptions have been made in some other studies of dielectric relaxation in proteins. [stI. [ssj It has 
been found[59j that there are only small quantitative differences between the dielectric relaxation associated with 
elliptical cavities compared to spherical ones. Clearly our results will be most relevant to globular proteins with a 
chromophore towards the centre. A more serious concern is that some chromophores can be located near the protein 
surface and so more exposed to the solvent. They may be better modelled by a spherical vacuum cavity at the 
planar interface between two different dielectric medium. Similar geometric considerations apply to transmembrane 
photosynthetic proteins and systems containing lipid membranes. The approach used here could be extended to such 
cases by considering continuum dielectric models with different geometries. (goI l6l| 

Point dipole approximation. The chromophore is treated as a point dipole and its spatial extent is neglected. More 
realistic treatments of spatially extended distributions inside solvent cavities, [13 do not lead to qualitatively different 
behaviour. 

Neglect of the nuclear dynamics of chromophore. Most electronic transitions will be associated with some structural 
change of the chromophore. Except, for the case of chromophores (such as retinal and PYP) which undergo confor- 
mational change upon photo-excitation, generally the reorganisation energy (and associated Stokes shifts) for modes 
with frequencies less than 1000 cm"-"^ are typically of order of lO's cm~^ and so much smaller than those associated 
with the solvent and protein. [6^. Furthermore, intra-molecular vibrations with substantial reorganisation energies 
have sufficiently high frequencies that they occur on timescales much faster than most of the experiments we consider 
and are not thermally excited at room temperature. 
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D. Overview of the paper 

In this paper we consider five distinct dielectric continuum models of the environment of a biological chromophore. 
For each model we derive an expression for the spectral density This allows us to explore how the relative 

importance of the dielectric relaxation of the solvent, bound water, and protein depends on the relevant length scales 
(the relative size of the chromophore, the protein and the thickness of the layer of bound water) and time scales 
(the dielectric relaxation times of the protein, bound water and the solvent) as is discussed in Section |Vl Many 
experimentally obtained spectral densities can be fitted to a sum of Lorentzians (see Table . For a protein that 
is large compared to the size of the binding pocket of the chromophore and the width of the bound water layer, our 
models predict a spectral density given by the sum of three Lorentzians which correspond to the dynamics of the 
protein, bound water and bulk water dynamics respectively. An essential feature is the separation of time scales, 
associated with the solvation coming from each of the three components of the environment. 

Specifically, to a good approximation, the spectral density is given by, 

_ apuj atuj a,uj 

'^^>' 1 + (c.r,)2 + 1 + {u;ny + 1 + (c.r,)2- 

The subscripts x = p,s,b refer to the protein, solvent and bound water, respectively. We show that when the dielectric 
relaxation of the different components of the environment are treated in a Debye approximation, that the relaxation 
times can be expressed as: 



2ep,, + 1 



(9) 



(10) 



Pi* 



n = TD,b (11) 

where €x,s, £x.iTTD,x are, respectively, the static dielectric constant, high frequency dielectric constant, and relaxation 
times of a Debye model for each medium, x ^ p, s, b. The high frequency dielectric constant is related to the refractive 
index n^, by e^^i = nl. 

We show that the reorganisation energies associated with each part of the environment are given by ax/T^, where 



Op _ SjAfiy {ep^s - £p,i) 

Tp neoa^ {2ep,s + l)i2ep,i + 1) 



Ctb 



3(A/i)2 (c-b\ (e^_, + 2ej_J(eb,,-efa,,) 



Tb 27reo53 \ b J el^^{2es^s + ^p,iY 
In particular, for typical systems the above three quantities can be of the same order of magnitude, i.e. 



as ab ap 

Ts Tb Tp 



(12) 



(13) 



(14) 



(15) 



Due to the large separation of time scales, the spectral density ([8]) will have peaks at approximately, uj = I/t^;. Hence, 
the peaks are of approximate magnitude, J(uj = I/tj,) ~ ax/Tx^ and can be of the same order of magnitude. This is 
because although each contribution involves different dielectric constants, they only have a limited range of values, 
and the ratios of the different dielectric constants that appear on the right hand side of the above expressions are all 
of order one. This is supported by experimental data (see Table where several relaxation times are observed which 
vary by several orders of magnitude, but whose relative contributions are of comparable magnitude. Hence, in many 
cases only a single component of the environment (protein, bound water, or bulk solvent) will be relevant to a given 
process. 
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The effect of the protein. The expression allows us how the ultra-fast solvation associated with the solvent is 
modified in the presence of a protein. Water is a highly polar medium with e^^s — 80, and e^^i ~ 4, and fast dielectric 
relaxation, td,s — 8 psec. Hence, we see that even if Sp^i is as large as 5 that the solvation time associated with the 
solvent is only be 50 per cent larger compared to the time of 0.3 psec which occurs in the absence of the protein. 
The effect of the protein on the strength of the coupling of the chromophorc to the solvent is more substantial. In 
(fT3|) . the coupling scales with the inverse cube of the power of the radius of the protein. Hence, if the diameter of 
the protein is four times the size of the chromophorc, the coupling as will be reduced by two orders of magnitude. 
The above results show that even a distant solvent can lead to ultra-fast solvation comparable to that found in the 
absence of the protein. This leads us to suggest that some studies which claim to have identified ultra-fast dielectric 
relaxation of proteins [63l. [6^ may in fact be detecting the fast response of the distant solvent. 

Order of magnitude estimates of a^- Typical values of parameters are a ~ 3A, b ^ 10 A, A/i ~ ID, and so the 
reorganisation energy axh/Tx ~ 10 meV ~ 100 cm^^. From Figure 2 and the Table in the Appendix we see that 
h/Tx is of the order of 10 meV, 1 meV, and 0.01 meV, for the solvent, bound water, and protein, respectively. Hence, 
the dimensionless couplings as ^ 1, ab ^ 10, and ap ~ 100. Hence, the only quantum dynamics that is likely to be 
coherent is that which occurs on timescales comparable to or faster than the relaxation of the bulk water, i.e., less 
than a picosecond. 

The outline of the paper is as follows. In Section [TT] we describe how the interaction between a chromophorc and 
its protein and solvent environment may be modelled by an independent boson model. We show how the interaction 
with the environment leads to decoherence of the electronic states of the chromophorc. In Section IIIII we propose 
a set of continuum dielectric models suitable for describing the environment around a chromophorc, and use them 
to obtain an expression for the spectral density in each case. In Section IIVI we consider particular limits of these 
spectral densities and obtain simple expressions for the contribution of each component of the environment (protein, 
bound water, and bulk water) to the total spectral density. In particular, we are able to obtain expressions that 
can be used to evaluate the relative importance of each component of the environment. We find that even when the 
chromophorc is completely surrounded by a protein it is possible that the ultra-fast solvation (on the psec timescale) 
is dominated by the bulk solvent surrounding the protein. In Section |V] we discuss methods for obtaining spectral 
densities from optical spectroscopy, and compare the predictions of our models to experimental data. In Section IVII 
we relate our results for the spectral density due to dielectric relaxation to what has been learnt from molecular 
dynamics simulations on specific protein systems. 



It can be shown [6| that the coupling of the electronic excitations in a chromophorc to its environment may be 
modelled by an independent boson model (65j , which has the Hamiltonian, 



We note that this is just the spin boson model with A = 0. Here the chromophorc is treated as a two level system with 
energy gap e between the ground and excited state. The first term describes the energy of the isolated chromophorc, 
described by Pauli sigma matrix cr^. The second term is the energy of the surrounding environment (protein and 
solvent), where the environment is modelled as a bath of harmonic oscillators [6j| . The final term describes the coupling 
of the state of the chromophorc az to the environment. In this case, the coupling [gl, H^l is due to electrostatic 
interactions between the chromophorc dipole and the "cage" of polarised solvent and protein molecules around it, as 
will be described in more detail below. The effect of this coupling on the quantum dynamics of the chromophorc is 
completely specified by the spectral density, defined by eq. (O |65| . 

If the two-level system is initially {t = 0) in a coherent superposition state l^t a\l > +b\2 >, which is not 
coupled to the environment. An ultra-fast laser pulse can create such a state. Then at time t, the TLS is described 
by a 2 X 2 reduced density matrix, p{t). It has matrix elements [68j 



II. QUANTUM DYNAMICS OF THE INDEPENDENT BOSON MODEL 




(16) 



Pii{t) 

P22{t) 
Pl2(i)-P2lW 




(17) 



where 9(t) is a phase shift given by 




(18) 
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and 

describes the decoherence due to interaction with the environment. 

The phase shift can be used to define a time dependent Stokes shift of the energy separation of the two levels. The 
instantaneous energy is found by taking the derivative of (|18p with respect to time, (Sa . [69j 

^(^t) = e - ^ = e-ER- r duo^ cos{ujt) (20) 
at Jq lu 

Hence, the spectral density can be determined by taking the Fourier transform of measurements of v(t), as discussed 
in Section IVl 

Depending on the relative size of the time t to the time scales defined by 1 /ujc (the relaxation time of the bath) 
and h/kgT, there are three different regimes of time dependence. For short times tUct < 1, 

r(t,r) = ^ (21) 

and so there is a Gaussian decay of coherence, on a time scale Tg, given by 

1 

Z2 ' 



d^J(^)coth( ) . (22) 



If in addition, ksT 3> hcuc, this reduces to 

n 

Tn 



= V^EnkBT (23) 



where Efi is the reorganisation energy given by Eqn. ([5]). 

For intermediate times, 1/ujc <t< Ti/ksT (the quantum regime 38]; it only exists if fc^T < TiWc), then 

r(i,r) »aln(Lj,0, (24) 

where a = J'{uj = 0), leading to a power law decay of coherence. 
For long times {t ^ Ti/kBT, the thermal regime) 

r{t, T) « 2akBTt/h = t/Td, (25) 

giving exponential decay of coherence. 

The crossover from Gaussian to exponential decay has been explored in quantum measurement theory in the 
context of continuous measurement and the quantum Zeno effect |47| . The decay of the off-diagonal part of the 
density matrix results from decoherence from the interaction of the TLS with the environment. Thus, we see the 
effect of the environment on the TLS is 



""-'Hit;" 'o'l^i'i 



Hence, the timescale Td can be interpreted as the "collapse" of the wavefunction of the TLS due continuous measure- 
ment of the state of the TLS by the environment. [H, [33, [13 We will see that using the spectral densities that we 
extract from experiment and from our continuum dielectric models that typically a > 1, and so at room temperature, 
Td < 10 fsec. 



III. THE SPECTRAL DENSITY FOR THE DIFFERENT CONTINUUM MODELS OF THE 

ENVIRONMENT 

In the simplest continuum model [STt picture of protein-pigment complexes, the chromophore can be treated as 
a point dipole inside a spherical dielectric [stI. | 67| . representing a globular protein, surrounded by a uniform polar 
solvent with complex dielectric constant es{uj) [66|. This can also apply to a chromophore-protein complex embedded 
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FIG. 3: The five continuum dielectric models considered for a chromophore-protein-solvent system. The chromophore is 
modelled as a point dipole. In Model 1, the chromophore is modelled sits at the centre of a cavity of radius a roughly of the 
van de Waals size of the chromophore, surrounded by a uniform polar solvent with complex dielectric constant ((.<;). In Model 
2, the chromophore is surrounded by an infinite protein, modelled as a uniform, continuous dielectric medium, with complex 
dielectric constant t^iiS)- In Model 3, the chromophore sits in a protein of radius h surrounded by the solvent. In Model 4, the 
chromophore sits in a cavity inside the dynamic protein, surrounded by solvent. In Model 5, the static protein is surrounded 
by a thin shell of bound water of radius c, surrounded by the bulk solvent. 



in a solid dielectric medium. In a previous work 'g'I, the spectral density was determined for a free chromophore in a 
solvent. However, many chromophores are inside proteins, which may have a significant effect on the coupling to the 
environment, if only in pushing back the solvent. 

Several continuum dielectric models for the protein environment have previously been proposed (see for example 
Figure 2 in ref. [STj and Figure 2 in ref. (70|.) These models often provide good qualitative and rough quantitative 
agreement [s^. We consider five distinct models, illustrated in figure[31 In every case, the central chromophore dipole 
polarises the surrounding cage of protein and solvent, which in turn produces an electric field inside the cavity, called 
the reaction field j66l . [67j . The interaction of this field with the central dipole is responsible for the interaction between 
the chromophore and its environment. The independent boson model can be obtained by writing the reaction field in 
terms of its normal modes, and quantising the coefficients in the standard second quantisation method Q, and the 
fluctuation dissipation theorem used to relate fluctuations in the reaction field to the appropriate spectral density 0] . 
A detailed derivation applicable to all models considered below, is given in the Appendix. 

Model 1. This describes a free chromophore with no surrounding protein. The molecule sits inside a spherical 
cavity of radius a, approximately the van de Waals radius of the molecule, inside a solvent with dielectric constant 
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es{uj)- The spectral density is|5£ 



Ji(i^) — 7^ rim— — . (27) 



where a is the radius of the cavity containing the chromophore and es(cij) is the complex dielectric function of the 
solvent and Cc the (static) dielectric constant of the cavity. A/x is the difference between the dipole moment of the 
chromophore in the ground and excited states. 

Model 2. This is somewhat analogous situation to Model 1, but in this case the chromophore is surrounded by an 
infinite, uniform protein with complex dielectric constant e.p{uj). This chromophore is again inside a cavity of radius 
a, which would be approximately the same radius as for Model 1. Such a model may be appropriate when the protein 
is very large. The spectral density has a similar form to Model 1, except involves the dielectric constant of the protein 
ep(w): 

J2{uj) ^ Tim / . 28 

Model 3. The chromophore is surrounded by a uniform dielectric sphere representing the protein. This protein 
sphere has radius b and (constant) dielectric constant ep. The spectral density is 

^ ' 27reo63 2e,(u;)+ep ^ ^ 

where es{uj) is the complex dielectric function of the solvent. 

Model 4. In a more detailed picture, (see the Figure the Appendix), Models 2 and 3 arc combined so that 
the chromophore sits inside a hollow cavity within the protein. The cavity has radius a (typically the size of the 
chromophore) and vacuum dielectric constant eo, while the protein again has radius b and is now described by 
complex dielectric constant ep{uj). Further detail may be added by treating the outer layer of the protein sphere as a 
separate, higher dielectric medium [STj representing the charged surface groups ^71|. In all cases, the chromophore is 
treated as a point dipole. The spectral density is given by 

' 27reoa3 2(ep - l)(e, - ep)a3 + (2ep + l)(2e, + ep)63 ^ ' 

where b is the radius of the protein containing the chromophore, a is the radius of the cavity containing the chromophore 
(usually the size of the chromophore), is the complex dielectric function of the solvent and ep is now the complex 
dielectric function of the protein. Note that the frequency dependence of the dielectric constants has been omitted 
for clarity. We see that for appropriate limits {{a/b) ^ 0, Cp = e^, Cp = 1, etc) Models 1 and 2 can be recovered, as 
expected. (To obtain Model 3, one would have to allow the centre cavity in Model 4 to have an arbitrary dielectric 
constant.) 

Model 5. The chromophore sits in a static protein (no cavity) and is surrounded by a thin shell of bound water 
with a different dielectric constant to the bulk solvent. To obtain the spectral density, we can use the results for 
Model 4, with ep — s- eb, Ec — ^p, a — > 6, and 6 — > c: 

U ^ = 1M1t (gfc + 2gp)(e. - + (£b - £p)(2£. + eb)c3 , . 

' 27reofe3 2(eb-ep)(e,-ef,)fe3 + (2eb + ep)(2e, + eb)c3 ^ ' 

Note that tp refers to a constant (typically high frequency) protein dielectric viz. Model 2. e;, is the complex dielectric 
of the bound water. 



A. Debye model for the frequency dependence of the dielectric constants 

To specify the dielectric constant of each component of the environment €x{i^) {x — s,p, b) we consider the Debye 
form of the dielectric constant [soj. 

e^iuj) = + ^""'^ . ^"^'^ (32) 

1 - lU}TD,x 

where td_x is the Debye relaxation time of the relevant component, ex,s = e(w = 0) is the static dielectric constant, 
and e^i = e(w oo) is the high frequency dielectric constant, within the range of physically relevant frequencies - 
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note that ex.i would always be 1 for sufficiently high frequencies [721 ■ For water at room temperature, Cg^s = 78.3, 
es^i = 4.21 and Tns = 8.2 ps [zl]. For comparison, for THF (tetrahydrofuran) the values are es^s = 8.08, e^^i = 2.18 
and td,s = 3 ps \74\ . 

In Model 4, the protein is treated as a complex frequency dependent dielectric with dielectric constant SpiuA, and 
as such is allowed to relax and respond to the chromophore. We will consider the case of a Debye dielectric [T^], but 
a more complicated model including multiple relaxation times is also possible [t^. Typical values for ep are between 
4-40 depending on which part of the protein is of interest [s^, [zl, |73| . Studies also suggest that charged groups on 
the surface of the protein can skew the average value of the protein dielectric, and it may be more appropriate to 
model such proteins as having an inner and an outer shell with two different dielectric constants. The high frequency 
constant, which is the only value used by many studies, is more difficult to determine but is generally assumed to be 
between 1.5 and 2.5 [2^. [TSj. Section [VII discusses determinations of £p{oj) from molecular dynamics simulations. 

The appropriate dielectric relaxation time of the protein may be different from the protein relaxation times, which 
can be of order miUiseconds [43], as there are processes (e.g., vibration of bonds) on the order of femtoseconds 
pl32. Table 3.13] which may contribute to the dielectric function. These may be missed in studies of the dielectric 
constant on the nanosecond timescale 
[80|). Molecular dynamics simulations 



79 and is perhaps unobservable for aqueous solutions of proteins (e.g., ref. 
77| suggest a protein dielectric relaxation time of 10 ns for a peptide, while 



vibrations may be of the order of 100 fs which may apply in certain situations. Other studies have found no single 
relaxation times, with relaxation processes occurring across the entire experimental range of 20 ps - 20 nsfsil. [s^. 

For Model 3, an appropriate value for the constant dielectric of the protein must be chosen, which will depend on 
the frequency range of physical interest. For example, for frequencies greater than l/rp, where Tp is approximately 
the protein dielectric relaxation time, the protein will be well approximated by its high frequency value. 

The hydration shell of hydrogen bonded water molecules surrounding the protein will have a dielectric constant 
different to that of the bulk solvent, and have a longer Debye relaxation time. 



IV. ANALYSIS OF MODEL 4 AND MODEL 5 
A. Model 4 - Dynamic protein and dynamic solvent 

The full spectral density (Model 4) describes the total coupling to the protein and solvent, but is too complex to be 
written analytically and because of the multiple timescales involved may be non-Ohmic, in certain frequency ranges. 
However, there are many cases where the use of one of the simpler descriptions (Models 1-3) would be preferable. For 
example if the frequency dependent dielectric constant of either the protein or solvent were not known, it would be 
useful to have a simple criteria to establish whether obtaining these parameters is worthwhile. If protein contributes 
negligible coupling beyond pushing the solvent back to a new distance, then specially obtaining its dielectric through 
experiment or simulation would be unnecessary. 

If either the solvent or the protein can be deemed unimportant, then simulations of the chromophores do not require 
their inclusion, saving valuable computational power. Conversely, if (for example) the solvent can be shown to have 
a significant effect on the chromophore, then treating the protein only will be insufficient and solvent effects must be 
included. Hence, we explore under what conditions the protein dynamics may be neglected. The following discussion 
assumes that the dielectric relaxation time of the protein is longer than that of the solvent, which is expected to be 
true in the vast majority of cases. 



1. Limit of a small chromophore surrounded by a large protein 

In the case of an "infinite" protein a/6 0, the full spectral density (Model 4) reduces to J2{^), the spectral 
density for Model 2. In this limit, the solvent can be neglected when compared to the protein. We might naively 
assume then that, since the ratio a/h appears cubed, when the protein is several times larger than the chromophore 
the above spectral density can be used and the coupling to the protein is far stronger than to the solvent. However, a 
closer examination of ([50]) shows that if ep{uj) ~ 1, the solvent will be more significant even for large values of h/a, and 
a more appropriate expression for the spectral density describes the chromophore in a cavity of radius 6 surrounded 
by solvent. Hence, we seek a more precise statement for when the protein can be ignored. The spectral density can 
be rewritten: 




12 



We can expand this expression in p- . Note that although this is a complex quantity, provided both the real 

and imaginary parts are small compared to unity we can use the Taylor series expression S2±^ « ^ _|_ (q, _ 5c)2;. 



Furthermore, both the real and imaginary parts of the prefactor of this expansion coefficient will be less than unity, 
and a sufficient condition for this expansion to be valid is that a/b be small. We find 



= J2(w) + Js{^) 



2es(w) + Ep 



9eo 



(2e, 



(34) 



where Js{^) represents the solvent contribution to the spectral density. Note that the dynamics of the protein only 
contribute to the first term, which conversely contains no reference to the solvent. The second term includes the 
solvent dynamics plus the high frequency limit of the protein dielectric constant as the only relevant protein property. 
Thus, we identify the first term with the protein contribution and the second term with the solvent, modified by the 
presence of the protein's high frequency dielectric. 



2. Relevance of protein vs. solvent 

By comparing the magnitudes of these two terms in p4p . we can establish the relative importance of the solvent 
and protein over different frequency ranges. We might expect that around cu — I/tj, and ui — l/r^ that the protein 
and solvent contributions should dominate respectively. Therefore, there should be a crossover point where the two 
contributions are roughly equal, somewhere in the range l/r^ <C w <SC 1 /t^ , which is what we look for now. 

As defined above, the cut-off frequency of uip — l/r^. Assuming a Debye dielectric for the protein, eq. can be 
approximated to first order in ujp/uj as: 

J2{uj) LJpUp / UJ , uj:s>ujp (35) 

Therefore, the tail end of the spectral densities given fall off as l/w (as compared to their linear rise for a; <C Wp). 
Noting that the reorganisation energy for this spectral density is Ep = apUip, we write the spectral density for uj :s> ujp 
as J2(w) = EpUjp/uj. (This is quite different to the spectral density for lo <^ ujp, which is J(w) = EpUj/ujp.) 

The second term in ([M)) is somewhat more difficult to evaluate. We will again Taylor expand in l/(Tpw). We note, 
however, that by extracting a factor of 1/a^ from both terms of (|34p the second term is proportional to (a/b)^, which 
we already assumed is "small" in the first Taylor expansion. Therefore, when we again Taylor expand around ujp/uj we 
need only work to zeroth order so that our total final approximation for Ji^uj) is to first order in two small expansion 
variables. This yields, 

Mc.) « ^!^J3(a;) (36) 

where J^iuj) is the spectral density for Model 3, with ep = ep^i, i.e., the chromophorc inside a constant dielectric (high 
frequency) protein with no cavity, surrounded by the solvent. The solvent contribution is therefore ohmic, with a 
dimensionless coupling constant, as given by (fT3|) . 

Therefore, equating the two spectral densities suggests a crossover between solvent and protein dominance at 
frequency Wco, 

^=.f^ (37) 

This ratio is always much larger than one for typical values of dielectric constants and relaxation times. 

For frequencies above this limit, provided we are in the regime 6 ^ a, we would expect that the protein dynamics 
are irrelevant for the system, and the dynamics of the chromophore is "slaved" to the solvent fluctuations. Similar 
effects have been observed in enzyme kinetics[83|. At low frequencies (w <C cJco), the protein dynamics dominate and 
the details of the solvent are mostly irrelevant. Hence, we expect that even when the chromophore is "shielded" from 
the solvent by the protein that the short time 1 psec) dynamics can still be dominated by the solvent. This raises 
questions about the recent assignment of the observed ultra-fast solvation to protein dvnamics. [63l [6^. [83| 



B. Model 5 - Bound water 



Our goal is to obtain analytic criteria which tell us when the bound water is relevant. If the dielectric contribution 
of the bound water dominates over that of the protein, we can use Model 5 to describe the system. Instead of 
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the chromophore pocket being treated as the cavity, now the entire protein is treated as a cavity of radius b with 
frequency- independent dielectric ep. This is surrounded by a shell of bound water with radius c (so the shell has width 
c — b) and dielectric €h{Lj) (with the subscript representing the bound water). We expect that the layer of bound water 
(typically flS\ about 4.5A) will be thin compared to the rest of the protein (6 ~ 20A), and so we are interested in the 
limit b c. 

We simplify the spectral density (l3T|) in the same way as Section IIV Al The full details may be found in the 



Appendix. Taylor expanding in ^-r^ yields 



(38) 



The first term represents the spectral density of a chromophore inside a cavity of radius b with dielectric constant 
ep surrounded by a bulk solvent, and so is the spectral density in the absence of the bound water, as described by 
Model 3. The second term, is proportional to the ratio ^j^, a measure of the thickness of the bound water, and can 
be identified with the contribution to the spectral density of the bound water: 

-^bwi^) 



2 

£6 



2e, 



(b - es \ (c-b) 



2e, 



£6 



Then, the bound water term can be expressed as 



(39) 

(40) 
(41) 

Using a Debye form for the bound water spectral density, gives Im[e6(tj)] = Ub s — i) ^^'^i -i ■ However, we must 

also include the |eb(ti;)p contribution to the frequency dependence. If we consider frequencies much less than the bulk 
solvent relaxation time I/tj, then we again find an ohmic spectral density for the bound water contribution, with 
dimensionless coupling given by Equation (jl4p . 
In comparison with the solvent contribution. 




Ob 



£6,1 



(42) 



We would typically expect Cs.s 3> es,i , ^b,s S> ct.i, and es^s 3> £6,^ and the protein static dielectric constant to be small 
compared to any static frequency, but perhaps comparable to the high frequency values of the dielectric constants of 
bulk or bound water. Therefore, 



£M-eM (gL + 2£s,J(2£s.. + £p) 



£fc.; 



4£ 



2£ 



£fc,i 

£fc,; 



which we expect to be of order one. Therefore, we would usually expect 

Off, C — b Tb 



(43) 



(44) 



' pe 



Hence, if t;, ^ Ts, as is observed |107l . ll08l |. then we would expect the bound water to be the dominant effect. Further, 
since the heights of the peaks are approximately given by their reorganisation energy {J{ujc) = aujc ~ Eji) we find 



•-b 



(45) 



where Eb, Eg are the reorganisation energies of the bound water and solvent respectively. 

The cross-over frequency between the bound water and bulk water contributions dominating the spectral density 
can be estimated by the condition JsitOco) = Jb{^co)- Assuming that the bulk and bound water timescales are 
sufficiently separated that at the cross-over point Jf,(w) is in the decaying tail and Je(w) is in the linear region, then 
the cross-over frequency is given by 
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V. SPECTRAL DENSITIES DETERMINED FROM ULTRA-FAST OPTICAL SPECTROSCOPY 



The spectral function J{uj) associated with optical transitions in chromophorcs can be extracted from ultra-fast laser 
spectroscopy [s^ IH, [s^ . Two widely used techniques for doing this are the dynamical Stokes shift and three-photon 
echo spectroscopy. 

Dynamical Stokes shift. The time dependence of the Stokes shift in the fluorescence spectrum, where is the 
maximum (or the first frequency moment) of the fluorescence spectrum at time t, can be normalised as 

Cit) = ^1^^ (47) 

such that C(0) = 1, and C(oo) = when the fluorescence maxima has reached its equilibrium value. Using Equation 
(PO)) this is related to the spectral density by 



dLu^^^coa{ujt) (48) 



where Ej^ is the total reorganisation energy given in eq. (O, which also equals the total Stokes shift associated with 
solvation. 

The function C (t) is sometimes referred to as the hydration correlation function and experimental results are often 
fltted to several decaying exponentials, 

C{t) = Ai exp(-i/n) + A2 exp(-i/T2) + A3 expi^t/rs) + ... (49) 

where Ai + A2 + . . . = 1. From (|48p , this corresponds to a spectral density of the form 

J w = ^ — + ^ — + . . . 50 

l-H(wTi)2 1 + (CJT2)2 ^ ^ 

The dimensionless couplings aj {j = 1,2,...) are related to the total reorganisation energy by 

— 2_J_I ~ 0.25A,- — (51) 
Trn. cm ^ psec 

Table HIl gives values of the fitting parameters {E^, Aj,Tj) determined by fast laser spectroscopy for a range of 
chromophorcs and different environments, both protein and solvent. We do not claim the list is exhaustive of all 
the published values, but is meant to be indicative (For example, see also [1^ IH, |63, HE Hlj ) We note the 
following general features. 

(i) The Stokes shift varies significantly between different environments, both solvent and protein. Generally, the 
presence of the protein reduces the total Stokes shift and the relative contribution of the ultra-fast component, which 
can be assigned to the solvent. The less exposed the chromophore is to the solvent the smaller is solvent contribution 
to the spectral density. This is also seen in measurements of the dynamic Stokes shift for a chromophore placed at 
three different sites in the Bl domain of protein G. (See Fig. 3C of Ref. lool ). Denaturing the protein tends to expose 
the chromophore to more solvent increase the total Stokes shift and increase the relative contribution of the ultra-fast 
component. 

(ii) The different decay times observed for a particular system can vary by as many as four orders of magnitude, 
ranging from lO's fsec to a nsec. 

(iii) The relative contributions of the ultra-fast (lOO's fsec) and slow (lO's psec) response are often of the same order 
of magnitude, consistent with equation (|15p . 

(iv) Even when the chromophorcs are inside the protein, the coupling of the chromophore to the solvent is large. 
For example, Prodan is in a hydrophobic pocket of HSA, well away from the surface, and yet as ~ 50. Even for the 
"buried" chromophorcs (Leu'' and Phe"^°) in GB1,[9^ the solvent contribution is AgEft. ~ 100 cm~^, Tg ^ 5 psec, and 
so as ~ 100. There are several proteins for which a very slow (~ lO's nsec) dynamic Stokes shift has been observed 
and has been assigned to dielectric relaxation of the protein itselffsil. [s^ l . 

(v) The values of the dimensionless couplings aj that we obtain from ([5T|) are comparable to the rough estimates 
we made in Section [I Dl 

Three pulse photon echo spectroscopy. This technique is analogous to stimulated spin echo measurements used in 
nuclear magnetic resonance to extract the phase relaxation time T2. For "long" times it can be shown[8^] that the 
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TABLE II: Solvation relaxation times for various chromophores in different protein and solvent environments. The values of 
relaxation times and their relative weights are determined by fitting the time dependence of the dynamic Stokes shift (|47p to 
the functional form (|49p . En is the reorganisation energy, given by (|6]), and equals the total Stokes shift. Note there is some 
variation in estimates of the reorganisation energy depending on whether one estimates it from the maxima in the absorption 
and emission spectra or from the first frequency moment of the spectra [3]. It should be noted that the time resolution is 
different in the various experiments. Some did not have access to femtosecond time scales and so we have left the relevant 
columns blank. SC is Subtilisin Carlsberg. HSA is Human serum albumin. SNase-WT is Staphylococcus nuclease in the wild 
type. SNase-KllOA is a specific mutant of Staphylococcus nuclease. HSA is in its native folded form in the buffer but denatures 
in concentrations of Gdn.HCl (guanidine hydrogen chloride) greater than about 5M and at a pH above 7. Trp is the amino acid 
Tryptophan, Acrylodan is 6-acryloyl-2-(dimethylamino)naphthalene, Prodan is 6-propionyl-2-(dimethylamino)naphthalene, AP 
is amino-pyridine, DCM is 4-(dicyanomethylene)-2-methyl-6-(p-dimethylaminostyryl) 4H-pyran, MPTS is 8-methoxypyrene- 
1,3,6-trisulfonate, and bis-ANS is l,l-bis(4-anilino)naphthalene-5,5' disulfonic acid. 



Chromophore 


Protein 


Solvent 


Ref. 


Er (cm"^) 


Al, Tl 


A2, T2 


^3, T3 


Trn 


none 


water 


[911 




0.65, 160 fsec 


0.35, 1.1 psec 




Trn 


none 


water 


[921 


2193 


0.55, 340 fsec 


0.45, 1.6 psec 




Trp 


SC 


buffer 


[931 


1440 


0.6, 800 fsec 


0.4, 38 psec 




Trp 


Monellin 


buffer 


[321 


960 


0.46, 1.3 psec 


0.54, 16 psec 




Trp 


SNase-WT 


buffer 


[941 


850 


0.46, 5 psec 


0.54, 153 psec 




Trp 


SNase-KllOA 


buffer 


[941 


876 


0.77, 3 psec 


0.23, 96 psec 




Trp 


HSA 


water, pH 7 


[951 


1156 


0.39, 5 psec 


0.61, 133 psec 




Trp 


HSA 


water, pH 9 


[95] 


1015 


0.3, 1.6 psec 


0.7, 46 psec 




Dansyl 


SC 


water 


[93] 


1180 


0.94, 1.5 psec 


0.06, 40 psec 




DCM 


HSA 


Tris buffer 


[82] 


515 




0.25, 600 psec 


0.75, 10 usee 


Prodan 


none 


buffer 


[96] 


2313 


0.47, 130 fsec 


0.53, 770 fsec 




Prodan 


HSA 


buffer 


[96] 


916 


0.19, 780 fsec 


0.56, 2.6 psec 


0.25, 32 psec 


Acrylodan 


HSA 


buffer 


[96] 


1680 


0.23, 710 fsec 


0.41, 3.7 psec 


0.36, 57 psec 


Acrylodan 


HSA 


0.2M Gdn.HCl 


[96] 




0.16, 280 fsec 


0.36, 5.4 psec 


0.48, 61 psec 


Acrylodan 


HSA 


0.6M Gdn.HCl 


[96] 




0.2, 120 fsec 


0.55, 2 psec 


0.25, 13.5 psec 


MPTS 


none 


buffer 


[62] 


2097 


0.8, 20 fsec 


0.2, 340 fsec 




MPTS 


Ab6C8 


buffer 


[62] 


1910 


0.85, 33 fsec 


0.1, 2 psec 


0.05, 67 psec 


bis-ANS 


GlnRS (native) 


water 


[33] 


750 




0.45, 170 psec 


0.55, 2.4 usee 


bis-ANS 


GlnRS (mohen) 


urea 


[33] 


500 




0.63, 60 psec 


0.37, 0.96 usee 


4-AP 


GlnRS (native) 


water 


[33] 


1330 




0.85, 40 psec 


0.15, 580 psec 


4-AP 


GlnRS (mohen) 


urea 


[33] 


700 




0.77, 50 psec 


0.23, 0.9 nsec 


Zn-porphyrin 


Cytochrome-c 


water 


[97] 


170 




0.4, 250 psec 


0.6, 1.5 nsec 
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time-dependent echo peak shift S{t), where t is the time delay between the second and third pulse, is related to the 
correlation function C{t), given in by 

S{t) = ^C{t) (52) 

where Tg is the decoherence timescale given in eq. (P^ . and associated with the "collapse of the wave function." 

The solvation dynamics of the fluorescein dye eosin bound to lysozyme in an aqueous solution was studied and 
compared to that for eosin in water without the protein[l5|. For both systems, ultra-fast solvation relaxation occurs 
in about 10 fsec and is assigned to bulk water. However, for the lysozyme-eosin complex a slower relaxation also 
occurred on the scale of 100 psec. This was assigned as predominantly due to water bound to the protein, mostly in 
the first hydration shell. This can be compared with dielectric dispersion measurements ^lOS] which suggest that there 
are two solvation relaxation times of 4 psec and 270 psec. A molecular dynamics simulation of lysozyme in an explicit 
solvent environment of 5345 water molecules found a single solvation relaxation time of 100 psec [98|. Jordanides et 
aLfm used the dynamic dielectric continuum model of Song and Chandler [6Q] to extract the spectral density, based 
on four different dielectric models. The full time dependence of the solvation was best described by a model which 
included the frequency dependence of dielectric constant of both the lysozyme and the water bound at the protein 
surface. These models for the lysozyme complex can be compared to our models if some simplifying assumptions are 
made. In particular, we need to treat the lysozyme protein as spherical with the eosin complex at its centre. Models 
I and II in Ref. [lH| then correspond to our Models 3 and 5, respectively. Models III and IV are approximately our 
Model 4, with the appropriate choice for the protein dielectric constant of the protein. 



VI. COMPARISON WITH SPECTRAL DENSITIES DETERMINED FROM MOLECULAR DYNAMICS 

SIMULATIONS 

For several specific proteins molecular dynamic simulations have been used to determine several quantities relevant 
to this work: the static dielectric constant of the protein, the frequency dependent dielectric constant, solvation 
dynamics, or the spectral density associate d with an optical transition in a chromophore or an electron transfer. 
B 113, S H H S [iM EqI Hoi [Toi] We hope that our work will stimulate further simulations of the spectral 
density for specific chromophores and proteins in an aqueous environment. To determine it one needs to calculate time 
correlations of the (reaction field) electric field at the location of the chromophore within the protein. Equivalently, 
the spectral density can be related to the fluctuations in the energy difference between the ground and excited states 
of the system j99|]. We now briefly review some of the results on specific proteins that are relevant to this paper. 



A. Tryptophan in monellin and water 



Molecular dynamics was used to calculate the time correlation function C{t) for trajectories of a few nanoseconds 

[l04| . For free Trp in bulk water C{t) was flt to a bi-exponential decay function with Ai = 0.86 ± 0.04, ri = 70 ± 10 
fs, and A2 = 0.14 ± 0.04, T2 — 0.7 ± 0.2 ps. For Trp-3 in the protein monellin C{t) was flt to a tri-exponential form 
with Ai = 0.66 ± 0.02, Ti = 70±10 fs; A2 = 0.22 ± 0.02, T2 = l.OiO.lps; A3 = 0.12 ± 0.01, = 23 ±2 psec. The total 
reorganisation energy was calculated as Er — 3200 cm^^. The two faster decays were assigned to the bulk water and 
the slowest component (T3 = 23 ± 2 ps) was assigned to protein dynamics including the motion of the chromophore 
within the protein. This assignment is consistent with the interpretation of NMR measurements but is different to 
that given in the associated experimental measurements 32| of the time dependent Stokes shift. The latter assigned 
the slower time scale {'^ 20 psec) to the dynamic exchange between water bound at the protein surface (the flrst 
hydration shell) with bulk water. 



B. Protein GBl in water 



The dynamic Stokes shift of a chromophore at the site of several different amino acid residues within the Bl 
domain of the protein G was measured. |90i The residues were replaced with an Aladin chromophore at sites that were 
"buried" (Leu'' and Phe^°), partially exposed (Trp^'^), and exposed (Ala^"'), to the solvent. The more exposed the 
site the larger the dynamic Stokes shift and the faster the relaxation. Motivated by these experiments, Golosov and 
Karplus performed molecular dynamic simulations for this 56 residue protein in a solvent of 6205 water molecules. [l05j | 
They calculated the time-dependent correlation function for the electrostatic interaction energy of the site residue 
with the rest of the system. This quantity should scale with the energy gap correlation function. For eleven different 
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sites the solvent coverage (defined as the ratio of the surface area of the residue that is accessible to the solvent to the 
surface area of the isolated residue) ranged from 5 to 45 per cent. The hydration correlation function was found to 
vary significantly between sites, but all contained components that could be assigned to ultra-fast decay (on the 100 
fsec and 1 psec timescales, due to the surrounding water) and much slower relaxation (on the lOO's psec timescale) 
that could be assigned to coupled hydration and protein conformational dynamics. However, there was no simple 
correlation between the slow relaxation time scale and the extent of the exposure of the site to the solvent, contrary 
to the correlation found by others. [lO^ 

C. Frequency dependent dielectric properties of an HIVl zinc finger peptide in water 

This peptide consisted of 18 amino acid residues and was simulated in a periodic box containing 2872 water 
molecules. [73] It was simulated for 13.1 nsec and exhibited a clear separation of time scales associated with dielectric 
relaxation of the different parts of the system. The water had a dielectric relaxation time of 7 psec, comparable to 
that for bulk water. Dielectric relaxation of the protein was dominated by a time scale of 4.3 nsec comparable to that 
found in simulations of other proteins and comparable to the time scale for rotation of the whole protein. The static 
dielectric constant of the peptide was estimated to be 15. 

D. Frequency dependent dielectric properties of ubiquitin in water 

Ubiquitin is a small globular protein compose d of 76 amino acids. It was simulated in a cubic box containing 
5523 water molecules for runs of 5 nsec duration. |lOl| Time dependent correlation functions (which are the fourier 
transform of the frequency dependent dielectric constant) could be fit to sums of two decaying exponentials with 
different weights and relaxation times. For the dielectric relaxation the three dominant timescales observed were 7 
psec, 2.6 nsec, and 1.9 nsec. These were associated with the bulk water, with rotation of the whole protein, and the 
bound water and side chains at the protein surface, respectively. Recently the same group extended the simulations 
to 20 nsec and also calculated the frequency dependent diele ctric constant of solutions of the proteins apo-calbindin 
Dqk and the C-terminal domain of phospholipase C-t1. [10^ 

VII. CONCLUSIONS 

The focus of this paper has been on the coupling of optical transitions in biological chroniophores to their envi- 
ronment. However, the approach and results presented here can be readily adapted to other transitions involving 
two quantum states which differ in the value of their electric dipole moment. Examples include intersystem crossing, 
non-radiative decay via a conical intersection, electron transfer and proton transfer. 

We hope our work will stimulate more work considering the following general claims, which this paper has elucidated. 

(i) A valuable approach to modelling quantum dynamics in specific biomolecular systems may be in terms of 
"minimal" models such as the spin-boson model where the system parameters and spectral density are extracted from 
experiment and/or quantum chemistry and molecular dynamics. 

(ii) Even when the active site of a biomolecule is shielded from bulk water, the latter can still have a significant 
effect on the quantum dynamics of the active site, especially if the time scale of interest is comparable to the solvation 
time scale associated with the bulk water. This can lead to solvent fluctuations dominating protein dynamics and 
function (ssj . 

(iii) The environment of the active site can be divided into three distinct components: the surrounding protein, 
water at the protein surface, and bulk water. The times scales associated with the dielectric relaxation of each 
component usually differs by several orders of magnitude and so each makes a unique contribution the coupling of the 
quantum dynamics of the active site to the environment. Furthermore, the relative importance of each component 
depends on how the time (or energy) scale of the quantum dynamics, compares to time scale of the solvation associated 
with each of the components of the environment. Table 1 compares the associated energy scales. 

(iv) The time scales associated with decoherence and the "collapse of the wave-function" in these biomolecular 
systems are experimentally accessible. Given the high tuneability of these systems they could potentially be used in 
fundamental studies concerning quantum measurement theory. 
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APPENDIX A: TIMESCALES 



TABLE HI: Timescales for various processes in biomolecules and solutions. The radiative lifetime of a chromophore is order of 
magnitudes longer than all other timescales, except perhaps protein dielectric relaxation. MD refers to results from molecular 
dynamics simulations. Of particular relevance to this work is the separation of timescales, Ts <^ rt <^ Tp (compare Fig. [2]). 



Process 


Timescale 


Ref. 


Radiative lifetime 


10 ns 


[42] 


Internal conversion 


lOfs 


[42] 


Bulk water dielectric relaxation 


8 ps 


[73] 


Protein dielectric relaxation (MD), td,p 


1-10 ns 


[77,20]J 


Ultrafast solvation in water 


lO's fs 


[31] 


Fast solvation in water, Ts 


lOO's fs 


[31] 


Solvation due to bound water, rj, 


5-50 ps 


[32] 


Solvation due to protein, Tp 


1-10 ns 


[33] 


Covalent bond vibrations 


10-100 fs 


[42] 


Elastic vibrations of globular regions 


1-10 ps 


[42] 


Rotation of surface sidechains 


10-100 ps 


[42] 


Reorientation of whole protein 


4-15 ns 


[101] 



APPENDIX B: SOLUTION FOR THE REACTION FIELD 



A change in the dipole moment of the chromophore leads to a re-organisation of the environment, which produces 
a reaction field acting back on the dipole. This is shown schematically in Figure U) 



The reaction field for these models can be obtained by a generalisation of the techniques in ref. [67[ as follows: The 
electric potential, 4>{x,y,z), satisfies Poisson's equation — — p/e, where e is the local dielectric constant of the 
medium and p is the local charge density. Away from the point dipole and surface boundaries, p = and we must 
solve Laplace's equation V^(/) — 0. At the dielectric boundaries (and in general), must be continuous and because 
there are no free charges, D = eE = —eWcj) is also continuous across the boundaries. 

Although the protein is spherically symmetric, because of the point-dipole, the system has only cylindrical symmetry. 
If the spherically symmetric electric potential in each concentric dielectric shell is given by (t>i{r), (j)2{r), . . we can 
expand (pi in terms of spherical harmonics [IJ : 



A.,„r" + %i ) F„(cos 



n=Q 



We consider explicitly the case where we have a cavity surrounded by a single dielectric shell inside a bulk solvent. 
The central cavity has radius a, dielectric ec and potential </)c(r, w), the shell has total radius b (thickness b — a), 
dielectric ep and potential (j)p{r,0) , and the bulk environment is described by dielectric and potential (/)e(r, w). 
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FIG. 4: Model 4 for the interaction between a chromophore and its environment. The chromophore is treated as a point dipole 
sitting in a cavity of radius a in the centre of a spherical, uniform protein which is treated as a homogeneous dielectric medium 
of radius b. The protein-pigment complex is surrounded by a solvent, typically water, which is again treated as a homogeneous 
dielectric medium though actual molecules are shown for clarity of explanation. The chromophore's dipole moment polarises its 
environment, which in turn produces an electric field, the "reaction field", which interacts with the chromophore. Fluctuations 
in the environment will translate to fluctuations in the chromophore's energy. 



We can then apply the boundary conditions: 



(r ^ cx)) (Bla) 

(/>^ = ^cos6l (Bib) 

((/'p)r=6 = (0e)r=6 (Blc) 

{(l)c)r=a = {(l)p)r=a (Bid) 



(Blf) 



The first condition is that the potential must go to zero at infinity. This means that all coefficients with positive 
powers of r must vanish, i.e., A^^n = for all n. 

The second condition is the field from a point dipole. As this is the only free charge in the cavity, this is the only 
source term (inverse power of r) that will contribute to the potential (j){r,9). Since Pi (cos 6*) — cos 9, only the n = 1 
term is involved. Therefore, -Bc,n=i = A* ^^id Bc^n^i = 0. (Nothing is said about Ac,„). 

The final terms describe the continuity of the potential and its derivative over the boundary. The first condition 
gives 

E [^r>,nb- + ^ P4cos^?) = Y: ^^"(cos0) (B2) 

n=0 ^ ^ n=0 

Because the spherical harmonics P„ are orthogonal, we can consider each term of this sum as being equal, so 

^P,nb" + ^-^ (B3) 
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In a similar way, the remaining boundary conditions can be applied to produce a set of linear equations on the ^4^ „ 
and and Bi ^. We have six boundary conditions and six variables (each, of course, a function of n) and so we are 
able to solve for all parameters. However, we are only interested in the field inside the cavity, and in particular the 
unknown part Ac,n- We find that all the Ac^n are zero except for n = 1. Thus, the potential due to the surface charges 
is given by (/)c,aurt — — XMJ'COS^ = — where we find 

. ^ ^ 2 (fp + 2£,)(ee - ep)a^ + (e^ - ej(2ee + £p)b^ 

a^2{ep-e,){e,-ep)a^ + [2ep + e,){2ee + ep)h-i ^ ' 

The actual electric field in the cavity due to the surface charges but not the dipole itself - the reaction field - is 
then R = Rz — —Vcj)^ ^^^^{x^y, z) — x/^z, which will be a constant throughout the cavity, parallel to the dipole, and 
proportional to the dipole moment /i. The spectral density describing coupling of changes in the chromophore state 
to this environment is related to the zero temperature fluctuations in the reaction field [6|: 

J(w) = (A/i)2Re j dt e"^^ (i?(t)i?(0))y^o (B5) 

This can be shown by writing the reaction field R{t) in terms of its normal modes' creation and annihilation operators. 
A/i is the change in chromophore dipole moment on the transition from the ground to excited states. The fiuctuations 
in the reaction field {R{t)R{Q)) are obtained Q from the fluctuation dissipation theorem, and are proportional to the 
imaginary part of xi'^) (|B4p above, yielding 

J{uj) = 2{^^lflm{x{uJ)) (B6) 



Note the use of zero temperature fluctuations is a mathematical derivation only, and provided the appropriate temper- 
ature parameters for the solvent and protein are used, the resulting spectral density is applicable to all temperatures. 



